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The major goals of this effort were; (1) to determine the conditions 
for the onset of Marangoni flows in molten tin, a model system for low 
Prandtl number molten semiconductor materials; (2) to determine whether the 
flov/s can be suppressed by a thin oxide layer; and (3) based on experimenta]| 
and mathematical analysis, to predict whether oscillatory flows will occur 
in the float zone silicon geometry in space, and if so, could it be 
suppressed by thin oxide or nitride films. 

Techniques were developed to analyze molten tin surfaces in a UHV 
system in a disk float zone geometry to minimize buoyancy flows. The 
critical Marangoni number for onset of oscillatory flows was determined to 
be greater than 4300 on atomically clean molten tin surfaces. Further work 
is necessary to attain the second and third goals. Analytical results 
indicate that the disk geometry is more stable than the cylindrical geometr\j 
and experiments may have to be conducted in space to attain the objectives. 
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Summary - Overall Project 


If Marangoni convection during float zone refining of semicon- 
ductor materials, such as Si, produces oscillatory flow, then it is 
nearly certain that space processing of such materials will not produce 
crystals of superior quality. The major goals of this project were (1) 
to determine the conditions for onset of oscillatory Marangoni flow in 
molten tin, a system which models the low Prandtl number molten semi- 
conductors and then, (2) to determine whether the oscillatory flow, as 
well as steady Marangoni flow, could be suppressed by extremely thin 
tin oxide films, on the order of 6 to 100 a thick; and (3) based on the 
tin results and mathematical modelling analysis, to predict whether or 
not oscillatory flow will occur in the float zone silicon geometry, and 
if so could it be suppressable by means of films such as oxides or 
nitrides at the nanometer level of thickness. 

The project involved; (1) Developing an experimental technique 
for float zone solidification which minimized bouyancy driven con- 
vection, and allowed in situ Auger analysis of the molten Sn surface. 
(2) Developing surface analytical techniques to measure tin oxide films 
down to monolayer thicknesses. (3) Mathematical modelling to guide the 
experimental research and ultimately address goal 3. 

The project was funded for 5.75 years at a level adequate to sup- 
port one full time experimental scientist (Max Noack), 10 to 15% of the 
F.l. (J. D. Verhoeven) and half support for a graduate student (C. C. 
Hsu) and around 10 to 15% of a consulting staff member. During the 
second and third year Dr. A. J. Bevolo consulted on problems concerning 



Auger analysis. During Che third year, work began with Dr. W. N. Gill 
and graduate student C. C. Hsu on mathematical modelling, and Dr. Gill 
was supported as a consultant during the final 1.75 years of the 
project, and Hsu was supported half time for two years. 

Following are a list of the accomplishments of the project. 

1) Techniques were developed to analyze molten tin surfaces in a 
UHV system. ^ It was discovered Chat the Auger CMA could be used to 
do ELS (electron loss spectroscopy) and that by so doing the minimum 
thickness at which the oxide film could be detected to be continuous 
was extended Co considerably lower values than with Auger analysis 
alone. 

2) A disk float zone (DFZ) geometry was developed^* which allowed 
bouyancy flows to be minimized (Bond number « 0.05) and allowed 
Auger/ELS analysis on flat molten surfaces. The technique was extended 
to a model system, NaNOg, and Marangoni flow was demonstrated.® 

3) Experiments were carried out in a UHV environment which showed 
that with the disk geometry, (Ma)gj^ for onset of oscillatory flow 

is greater than 4300 on atomically clean molten Sn surfaces (unless the 
frequency of oscillation is greater than around 10 Hz).® 

4) Both a finite difference analysis and an approximate boundary 
layer analysis have been carried out for the DFZ geometry.®”® The 
finite difference analysis shows that for conditions of the Ma » 4300 
experiment, surface velocities of 25 cm/s are expected; quite good 
agreement was found between the boundary layer analysis and the finite 
difference analysis. 
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5) The boundary layer analysis was extended to other geome- 
tries®* ^ and it was found that for cylindrical geometries multiple 
solutions occur at critical values of the parameter 2(d/i)® (Ma/Pr), 
which suggests flow instabilities; whereas for the disk geometry, no 
multiple solutions were found. This result suggests that the disk 
geometry may be considerably more stable to oscillatory flow than the 
cylindrical geometry characteristic of float zone configurations. 
Publications Resulting from this Project 

1. J. D. Verhoeven, M. A. Noack and A. J. Bevolo, Suppression of 
Marangoni Convection with Oxide Films, Processing in the Reduced 
Gravity Environment of Space, Materials Research Society Symp. 
Proceedings, Vol. 9, 513-522, Ed. by Guy Rindone, North Holland, 
N.Y. (1982). 

2. A. J. Bevolo, J. D. Verhoeven, and M. Noack, Electron Loss Study of 
Native Oxide of Tin, J. Vac. Sci. Tech. ^(4), 943 (1982). 

3. A. J. Bevolo, J. D. Verhoeven and M. Noack, A Low Energy Electron 
Loss and Auger Study of the Oxidation of Liquid and Solid Tin, 
Surface Science 134 , 499 (1983). 

4. W. N. Gill, C. C. Hsu, J. D. Verhoeven and M. Noack, Thermocapil- 
lary and Buoyance Driven Flows in Liquid Metal Disks. Proceedings 
of PACHEC '83, Vol. I, 138-44 (1983). 

5. C. C. Hsu, W. N. Gill, M. A. Noack and J. D. Verhoeven, Marangoni 
Flow Velocity in a Thin Disk of Molten NaNO^, Chem. Engr. Comm. 
167-376 (1983). 

6. W. N. Gill, Floating Zone Crystallization, Symposium on Separation 
Technology, 206-210, Taipal (1983). 


w. N. Gill, Surface Tension Driven Boundary Layer Flow In Floating 
Zones, Chera, Engr. Coram. 57 (1983). 

W. N. Gill, N. D. Kazarinoff, C. C. Hsu, M. A. Noack and J. D. 
Verhoeven, Thermocapillary Driven Convection in Supported and 
Floating Zone Crystallization, COSPAR, Advances in Space Research, 
Pergamon Press, (1984 or 85). 

W. N. Gill, N. D. Kazarinoff and J. D. Verhoeven, Convective Dif- 
fusion in Zone Refining of low Prandtl Number Liquid Metals and 
Semi-conductors, Symposium on Chemical and Physical Processing of 
Integrated Circuits, ACS Symposium Series (1984 or 85). 


Summary - Incompletcd Work 


Unfortunately^ the project never accomplished Its major goals of 
establishing (Ma)gj^ values for onset of oscillatory flow in molten 
Sn and of determining the effectiveness of tin oxide on suppressing 
both oscillatory and steady Marangoni flow. It is felt that with an 
additional year or two of funding the experiments could have been 
extended to Ma values greater than 10, 000, where Schwabe has observed 
oscillatory flow in cylinder geometries of a Pr •• 9 fluid. Also, the 
effect of oxides on macrosegregation, and therefore on Ma flow, could 
have been evaluated, using radioisotope techniques; and further theo- 
retical work could have been performed on the influence of adsorbed 
oxygen and oxides on Marangoni flows in liquid metals and semicon- 
ductors . 

The contract terminated before a stability analysis could be per- 
formed on the nonlinear multiple solutions which were obtained for the 
base flows. It has been shown clearly in the paper listed as Reference 
9 above (attached in full, starting on page 18) that the pressure dis- 
tribution in the zone is intimately related to both nonlinear inertial 
effects in Marangoni flows and to the surface deflection, and this is 
conjectured to be important in hydrodynamic and thermal stability. The 
analysis of surface deflection for a supported cavity or slot, which 
appears in this paper, applies in large measure to the disk geometry 
also. It shows for the first time that the surface deflection, the 
pressure distribution, the parameter Re A/ c, and the parameter A^Re 
(which reflects the importance of nonlinear flow effects) are inti- 
mately related. It is conjectured that this interaction has an 


important influence on flow stability in the zone and may explain why 
existing stability analyses predict critical parameter values an order 
of magnitude too low. 


Suaunary of Work done October 1, 1983 to June 30 > 198A 

Wc were Informed ac the beginning of this year, October 1983| that 
the project would be terminated after nine months^ In June 198A> With 
this limitation In mind, and given the goals of the project which are 
listed on page 1, we decided to concentrate the experimental work In 
the following three areas; 

1) Attempt to determine the onset of oscillatory Marangonl flow in 
the DFZ apparatus, using miniature thermocouples to detect oscillatory 
flow. 

2) Attempt to detect changes In solid-liquid interface shapes in 
the DFZ samples due to presence of Marangonl flow. 

3) Attempt to determine the effect of Marangonl flow upon macro- 
segregation during plain front solidification experiments In the DFZ 
samples. 

In summary, the work progressed as follows: The experiments In 

the second area were unsuccessful because of an unanticipated effect 
which Is described In the following report. The experimental work In 
the first area was successful, but showed that we were limited to a 
maximum Ma number of A300 with our present heater, at which value no 
detectable oscillations occurred. Because development of a new heater 
would tak“ too long a decision was then required on whether to push the 
present heater, beyond design limitations and risk burnout, or to start 
the third area of experimental work. The decision was made to start on 
the third experiment and it was only partially completed by project 


termination on June 30 
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The raojor reeultB of the final nine montha work involve expori- 
mencal results from the first oroo and theoretical modelling work, both 
of which arc presented in our July 1984 COSPAR paper. In addition, the 
theoretical modelling work (including some work done on non NASA 
funding since the June termination) has been preftired for an ACS review 
paper. Hence, the report which follows includes essentially 3 


sections: 

I. July 1984 COSPAR paper — 9 

II. ACS Review paper 18 


III. Details of Experimental studies during this 
reporting period, October 1, 1983 to June 30, 1984 


Experiment I, Detection of Oscillatory Flow ------- 4 x 

Experiment 2, Interface Shape Evaluation -- - --50 

Experiment 3, Macrosegregafton Study --53 


I. Paper presented at the July 1984 COSPAR Conference in Graz, 


Austria and to be published by Pergamon Press (1984 or 85) • 
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J.D, Varhoevont 
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ABSTRACT 

Narangoni flow wai induced in a chin diak of uolcen tin with an atonically claan aurfaca by 
heating ic from below in on ultra high vacuum apparatua. The radial diatribucion of taapero-* 
cure waa beaaured and no tenperacura oaeillatlona were obaerved. Strong aechanletl diacur<- 
boncea introduced at a Harangoni number of about 4300 died out quickly with time indicating 
that the Byacea waa arable, Katheoacical Dodola Indicate that atrong flowa In a two cell 
atruccure exiat under the condiciona atudied, and that Che velocity profile in the call near 
the aurface la of the boundary layer type which variea vapidly with depth, 

INTRODUCTION 

The 'iA'p-irtonce of the developsenc of high quality alngle cryacalllne illicon in the evolution 
of tiu) :'>,rrent "high-tech" aolid atato Induacry ia clear, A aignificant llsiicaCion to the 
qv.eri,„' Jf alngle cryatal silicon la the presence of growth seriacions, which result from 
)rU>illatory convective flows present in the molten Si at the solid-liquid interface during 
irowth of the Si crystals from the melt. Because of natural convection, it is extremely 
difficult to eliminate cheao oscillatory flows at 1-g, but In space one expects no 
oscillatory flows from the buoyancy force. Hence, the preparation of Si single crystala free 
of strlaclon defects it a potential application for space prouesaing. Although ic. Is not at 
all clear that practical economics would favor such an operstioa, it is very difficult to 
predict these economics now for the long range future with the advent of such technologies as 
a space station, Therefore, it seems important that the potential of space processing for 
the preparation of scriatlon free semiconductor materials be established. 

Several experiments /I, 2, 3/ have shown that oscillatory Harangoni convection occurs in high 
Prandcl number liquids. At present there ie little agreement concerning the proper dimen- 
sionless parameter to be uecd to characcorlzo the critical physical conditions for the onect 
of such oscillatory flow, and there are no experiments known to us which have deteonatraced 
oscillatory Ktrongonl convection in low Prondtl number liquide such ae molten metola end 
semiconductors. If Harangoni convection during floac zone refining of semiconductor 
macerials euch as Si produces an ooeillatory flow, then it is nearly certain thnt spacn 
proccaslng of such materials will not produce the hoped for scrlation frue crystals of 
superior quality. Hence, the question of the conditions required for the onset of Harangoni 
flow in low Prondtl number liquids has importont practical implicacions, aa well ae interest- 
ing theoretical Impllcaclone, This paper preeents results of some initial experimental and 
theoretical atudiet an this problem. 

The zone studied both experimentally end theoretically is a circular disk of pure tin hcatod 
from below by a tontalum heater and cooled around the periphery of the disk by an oil bath 
whose temperature was controlled. The disks, which were approximately 1 to 1.2 mm in height 
■nd IS mm in radius, vested on the tantalum heater which varied from 1.7S to 3 mm in radius. 
By varying the power to the heater and controlling the temperature of the oil bath, molten 
zones of tin with radii from about 4 to 6.3 sn were produced. 

Differences between the high temperoture at the center of the disk and the solid-liquid 
interface at the periphery of che molten zone ranged from about SO to 140*C, and average 
temperature gradients across che disk varied from approximately 113 Co 230*C/cm, 

Hachematical Hodcl and Analysis of Surface Tension Driven Flow in Liauid Hetal Disks 


Two approaches have been taken here to model the surface tension driven flow in a liquid 
metal disk. In the first we have used a finite difference analysis of che Navior-Stokea and 


•ntrgy aquatlona and boundary condltlona which cloatly rtpraaanC tha raal phyalcal ayateis 
under ecnaid»raclon> In the aecond, wore alwpllfled approach, we have neglected "end 
elfecta" and have aaeuiaed the liquid region la unbounded In the radial direction which 
enabloa one to reduce the Navler-Stokee equatlona to an ordinary differential equation. 

Theeo approachea arc coaplccientary in that both of thca retain the eeecntlal nonlinear 
fKaturas of the problem, and the almpllflad approach auggeaca waya that the more detailed 
finite difference calculatlone can be correlated and Interpreted. The almpllflcd wodel alao 
enablea one to aearch for multiple ateady atate aolutlona of the problem and thereby get 
aonie Indication of the affect of the phyalcal and geometrical paromatera on the hydrodynamic 
atablllty of the Karangonl flowa, which occur In the matorlalo and conflguratl^ina of 
Intereat In cryatal growth. 

The equatlona of motion and energy for a dlak heated uniformly from below with a flux q 
over the radlua a, and syemotrlcally with reapect to ita axla are 
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The coordinate r la meeaurnd from tho axla of Che disk and x la measured from the upper 
surface down into the fluid, the total depch of fluid being d. 


The boundary conditions are 
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where i ig the radlua of Che molten tone end la Che melting point. 

Flnito difference aolutlona to equatlona (^-12) were obtained for Frandtl number, Pr - v/o, 
equal 0.014, the Pr of tin, and for Pr ~ 0.02, 9, the latter being the Pr of NoN 03 . To check 
the accuracy of the calculaclcns by being able to compare them with experimental data some 
solutions with Pr " 9 were obtained for a cavity and these reaults compare very favorably 
with Che surface velocltiea measured by Schwabs and Schsemann /!/ for NaK 03 . The agreement 
between our calculations and the experiments la especially good above a 4T of 35 ’C In which 
case the entire cavity was liquid. 

Figure (1) ehowB a comparison between the temperature dletrlbutiona predicted by aasumlng 
only diffusion In the melt <no convection), and a conatant heat flux over 0 s r 1, and 
chose obtained from the finite difference enalysla of equations (1-4), which Includes both 
conduction and convecdon. One aens that tha oure conduction solutions differ rather little 
from the temperature distribution produced by Marangoni convection until the value of Ha 
exceeds 2500. Thus, because of the high thermal conductivity of liquid metala, one can 
obtain a reasonably good representation of the temperature distribution by neglecting con- 
vection completely unlese Ms la large. This suggests chat considerable information about 
Marangoni flow may be obtained by using a simpllflad surface temperature distribution in 
solving the equations of motion one Is Interested primarily In low Frandtl number fluids 
Buch ae liquid ssetala or aeml-conductors. 

If view of the preceding dlscusolon, we decided to simplify markedly the problem aa follows: 
In sedition to assuming chat the free surfaces are flat, wp. assumed that surface heating Is 
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OF POOR QJAU«V 


•cconpliihad with « he«c flux which will produce • parabolic temperature dlatrlbutlon on the 
■urface of the dlak. Therefore the temperature dlitrlbuclon on the dlak aurface la given by 

“ ’'’hot " ('^UOT ■ ’’^cold) ff) 

One can ahow that as Pr 0| a constant heat flux on Che surface will produce precisely the 
distribution given by equation (13) t In equation (13) • TyOT corresponds to the caaperaturs 
at r " 0 and Tqoj^q represents chat at r " t. That la, 1 la simply Che distance over which 
one los/tsures the temperature difference dTi Equation (13) is Che pure conduction solution 
plotted In Figure 1. 




Fig, 1. Comparison of pure conduction Fig. 2. Surface velocity distribution for 

with a/1 1 and convection with a/1 ■ disk with 1/d ■ 6. Comparison of simplified 

0,875 aoluclons for disks with 1/d 6 model with finite difference solution for 

end various Karangoni numbers, various Ma, 


Ve further assume that "no end effects" exist. This means that the flow la not obstructed by 
a solid surface «c r » 1. With these assumptions , one can obtain exact aoluclons of equa- 
tions (1-3, 5-7) as follows. Define a stream function as usual by 


-ru-|| 

(14) 

and let ij' cake the form 


^ (o ^ 

(15) 

where n • x/d and 1^ a function of r only. One can show readily 

u ■ -fd dvg/dr and thus equations (1-3) are reduced to 

choc V ■ Vj^f and 
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(16) 

with 

f(0) - f(l) - f'(l) - f”(0) + 1 ■ 0 

(17) 

for a supported disk and 


f(0) - f(l) - f"(l) - f"(0) + 1 - 0 

(18) 

for a floating disk If the velocity scale, v^, Is 


2vfdV Ka r 
''r d llj Pr 1 

(19) 

and If d Is interpreted os the thickness (or height) of the supported disk and Che half- 
thlckneas of the unsupported disk, and v Is Che kinematic viscosity. The coordinate x Is 
mens’seec' from the free surface Into the fluid, and gj^ ■' g,. - 0, 


Tf,e quantities on the lefthand side of equation (16) ace functions of n only, and chose on 
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the right arc only of r. Therefore B !■ e constant which can be determined with the aid of 
tho boundary conditions given In equations U7} and (IB). Note that solutions of equations 
(16) to (IB) depend only on the single parameter 2(d/l)^Ka/Fr and computations were carried 
out for supported and floating disks up to values of the paranater equal to 90,000. 

Figure 2 shows some comparisons between the finite difference calculations of surface 
velocities and chose predicted for the supported disk by equation (16) with the boundary 
conditions (17). It Is seen that the agreement la excellent over BOX of the disk and in 
this case, for the finite difference computations , the disk was hected uniformly only over 
67. 6X of Its radius as Indicated by the value of a/1 at the bottom of Figure 2. Thus, the 
end effects are due to both Che solld*llquld Interface and the fact that the heater does not 
cover the entire liquid region as Is assuiud In the simplified suidel. Also, It Is clear 
that the agreement Is better at lower values of Ma, where eo;sductlon predominates and the 
simplified model results are essentially indistinguishable from the finite difference solu- 
tion. Note that the finite difference solutions show a double peak structure at larger 
Values of Ms, as shown by Che solid lines In Figure 2. 

The surface velocity distribution obtained from this elmpltfted model Is given by Figure 3 In 
which V[^ Is given by equation (19). Clearly, vg/vj^ - f'(0), Is a function only of 
(d/B)’Ka/Pr as Is the local velocity field when V]^ le used as the velocity scale. It is 
noteworthy that the vg/vg ratio for both floating and supported disks Is essentially the same 
except at small values of 2(d/t)^Ha/Pr, when d Is interpreted ac Indicated. 
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Fig. 3. DlmenslonlesB surfsce, velocity. 
Vg/vj^ as a function of 2(d/t)^Hs/Pr for 
supported and floating disks. 



Fig. 4. Dimensionless steady state radial 
velocity profiles for surface tension driven 
flow In disk supported on bottom. 


Radial velocity profiles as functions of distance from Che surface of the disk are shown for 
the simplified model In Figures 4 and 5 for supported and floating disks, respectively. 

The value of the parameter 2(d/i.)^Ka/Pr Is varied by two orders of magnitude to demonstrate 
how the shape of the velocity field changes with changes In either geometry, d/1, or heat 
flux, Mb. One sees chat the depth, over which Che velocity changes from the value at the 
surface to tero decreases from about 26X - 30Z of the disk height to about 7Z - lOX as Che 
characteristic parameter Increases from 200 to 26,000. Furthermore, the velocity in most of 
the core of the floating disk Is sloiosc constant ac higher heating rates and as Che depth to 
radius of the molten zone Increases. 




Fig, 6. Dimensionless steady state radial Fig. 6. Steady state scluclonp of Navler- 
veloclcy profiles for surface tension driven Stokes equaclonn for surface tension driven 
flow in floating dlak. flow in floating (1) and supported (2) disks. 
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Figure 6 llluscraces how the dioeitslonlese preeaure gradient paraa«cer> it (which In a con> 
■tent for a given value of 2(d/t}^Ma/Fr) varlea with the characterlatlc paranetert 
2(d/l)^Ka/Pr. One sees that B varies very slowly beyond values of 5000 for the characteris- 
tic paraoecer and at values of 90,000 B equals Itl x lO**^ foe the supported disk and -1,22 x 
10*^ for the floating disk. Over the range of 0 to 90,000, no evidence of multlplo steady 
state solutions to equations (16-18) waa found for either the supported or floating dlska. 

The radial velocity profiles all exhibit only two cells, one radially outward and one 
radially inward, for all values of 2(d/l)^Ka/Pr, This rather stable behavior is in 
dramatic contrast to that which we have found for the cavity and the cylinder, both of which 
evidence multiple aCeady-scate solutions and are less stable configurations from a fluid 
mechanical viewpoint, 

Exnerlmental Results 

Molten tin was chosen for study because its low Frandtl number, O.OIA, and its tendency to 
oxidize are typical of molten metals and semiconductors; and because Its low vapor pressure 
makes it amenable for study In a vacuum. It Is generally ogreed that Che strong tendency of 
molten metals and iiemiconductors to form oxide films on their surfaces will have s pronounced 
effect on Marangonl convection. Therefore an experimental design was daveloped which 
employed ultra high vacuum (UHV) techniques incorporating both an argon Ion sputter gun to 
clean the molten tin surface and on Auger analyser to examine the molten tin surface ^ situ 
and verify chat it is free from oxides. The system utilizes a disk float zone (DFZ) 
geometry which is shown schematically in Fig, 7, A tin disk of diameter 2R and height d Is 
cooled around Its periphery by a copper cylinder whose temperature is controlled by flowing 
hot oil supplied from an external constant temperature bath. The tantalum heater at the 
bottom center of the disk produces a molten zone of diameter 21 in the disk. The entire 
system la within a UHV chamber capable of achieving a vacuum of 2 x 10“^® Terr with the 
system at temperature. 


CYUiNDRICAL AUGER ANALYZER 



Fig, 7, Schematic representation of the disk float zone configuration. 

Two types of experiments were performed on the tin disks. In the first type, temperature 
traverses were made using a miniature thermocouple, which in some cases was W-Re with 50 pm 
wires, the bead of which was Inserted directly in Che molten metal and in ocher cases 
Involved a 25 mm OD thermocouple which was chromel-alumcl (type K) and was contained within 
a stainless steel tube placed with pure chromium, which was slightly oxidized to avoid 
reaction with the molten tin. 


Accurate experimental temperature distributions were obtained from these traverses, and two 
typical ones are shown in Figure 8. In Figure 8, we also Include temperature distributions 
calculated from solutions for the profile of temperature averaged over the height of the 
disk, T, which are given by 
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onmimL p .; . ;, 

OF POOR QUALITY 


Equacloni (20) »nd (21) caiUBM no r«dl«clon loia froa the tin to the vecutm lurroundlng it 
beceute thie can be ihoim to be email. 

One leee in Figure 8 that the temperature dlatclbutlon given by equatlona (20) - (22), In 
which the thermal conductivity of the liquid phaae, k|;^, la 0.32, J/cm*C, la conalatently 
higher than the experinental reeulta. Note that the heac flow la pTedoninantly radial in the 
dlak, and we Interpret the dlfferencea obtained In the esperimanta between the centerline 
temperaturea and the melting point being emeller chan choae calculated by equations (20 - 22) 
CO be strong evidence chat Karangonl flow exists In the liquid melts studied. However, the 
temperatures measured were very stable aa predicted by the elnpllfled models, and there le no 
evidence of oscillseory flow. Because very thin disks ware used, the maximum Bond numbers 
for our experiments was about 0.05. Therefore, natural convection waa negligibly small as 
confirmed by a finite difference simulation of run 56. The second type of experlisenc was 
conducted to determine If we could Induce thermal oeelllatlons in the system by perturbing It 
strongly after we had created the largest temperature difference, disk size and Harangonl 
number that we could on earth without haying the molten zone fall out. This experiment^ was 
run as follows. After achieving UHV conditions, the oil sink temperature was adjusted to 
166*0 as measured by t.c. f2, and a molten tone with a diameter around 8 mm was achieved. 

The Argon Ion sputter gun then was turned on and Che molten surface was cleaned by moving 
the beam over It. Both Auger and Low Energy Electron Loss (LEEL) signals were used to be 
sure that the surface was free of oxide. The thermocouple within the molten zone was 
monitored on a atrip chart recorder to detect possible oscillations resulting from oscilla- 
tory Harangonl flow. This procedure was repeated with Increased power and the maximum 
lulteu tin diameter achieved In this experiment was 11.7 ss>. The thermocouple tip was 
positioned at a radius of 3 nn. 

Temperature oscillations were not detected in this experiment which indicates that oscilla- 
tory Harangonl flow probably was absent. There are four important questions regarding this 
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Fig. 9. (a) Recorder trace of thermocouple FI 

with 21 - 11.7 mm and T- - 166"C. 

(b) Response of thermocouple output 
to a mechanical tap on the system. 



conclusion which will now be addressed. 


t 



Fig. 8. Comparison of pure conduction 
solution with experliaental data. 


Minimum amplitude detection limit . Figure 9 presents a recorder trace from the experiment 
which demonstrates the amplitude sensitivity of the experimental system. We estimate chat 
oscillations of amplitude greater than or equal to O.A uV would be detectable if present in 
the signal. . For the type K thermocouple used here, this would correspond to 0.01°C. 

The ability of the system to detect temperature oeelllatlons was qualltacively confirmed by 
giving the system a sharp mechanical tap with the result shown in Fig. 9b. Furthermore, the 
attenuation of the oscillations produced by tapping the system Indicates its stability to 
finite amplitude perturbatlona. 

Frequency response . The expected frequency of Harangonl Induced oaclllstory flow in low 
Prandtl number systems in the disk geometry is unknown. If the frequency were high, the 













omGm.L 

OF POOR » V 


United cine reiponae of our recording eyecen would reduce Che uiplltude of elie ncaoured 
oeclil«clon> The cine responae of our recorder was evaluated by neaauring Che decreaee In 
anpllcude of a alne wave aupplled by a frequency lenerator aa a funccion of frequency* 
uclllzlng Che aane aeale aa In Che experlnenc. In addition co llnlcatlona rcaulClng fron 
Che recorder, one nuac alao conaider che loaa of anpllcude due co the flnlce Cine reaponae 
of the thermocouple Icaelf. The tine conatant of a 0.2S vn dlaauiter thermocouple, t, haa 
been ncasured /&/ In molcen So-Pb alloya to be 35 na. The fracclonal drop In amplitude aa 
frequency, f, Increaaea can be taken /5/ aa, 1/(1 + (2af)^)*'^, The total fractional lota 
of amplitude was taken aa the product of that due to tha recorder and the thermocouple, and 
che reaulti ahov chat che minimum detectable aapllcudea of the temperature oaclllatlona are 
0.01, 0.013, 0.062, 0.34 and 3.15*C for frequanclaa of 0.1, 1, 5, 10 and 20 Bz. Thua, 
camperafure oaclllatlona of anplltudea greater chan 0.5'C could have been praoant In the tin 
dlaka and would not have been detected if their frequenclea were above around 10 Bz. 

Cleanllneaa of molten tin aurfaee . In prevloua work /6/, we have ahown Chat one may utilize 
LEEL apeccroBcopy with the cyllndrlcel mirror Auger analyeer to detect tin oxide layera on 
cnlten tin of less chan a monolayer. By utilizing this technique, coupled with Auger 
apeccroscopy co evaluate other surface contaminanta, the following picture emerged. 

The original tin surface was carefully eleccropollahed. Before the initial melting at a 
pressure of leas chan 10"® Torr, the surface wae contaminated by a very thin C layer on top 
of tin oxide. When this cnntamlnatlon was removed by sputtering before melting, it waa 
found that upon molting the molten surface became contaminated with both oxygen and sulfur. 
These contaminates could be removed by moving che argon epucter beam over the molten 
surface, and the surface would remain clean over 12 hour periods. 

However, If che power were increased to melt out to a larger zone radlua, che aurfaee became 
contaminated again. Obviously, oxygen and aulfui* Impurities In the solid tin contaminated 
the surface upon melting. We were never able to clean the surface completely out to Che 
full molten zone radlua, perhaps due co surface diffusion. By moving the sample under the 
Auger analyser, It was detemined that at the maxistua molten zone diameter of 11.7 mm, the 
central 10.6 mm waa atomically clean, but an outer annulus of 0.65 mm width was covered with 
a few monolayers of contaminant. It seems doubtful that this small amount of contamination 
would affect the Marongonl flow, particularly since the maximum temperature gradients occur 
in the atomically clean part of the zone. 

Evaluation of Marancortl and Reynolds numbers . Schwabe, et al. /!/ have presented careful 
experimental studies which detect the onset of oscillatory flow In NoN 03 (Pr ■■9) at a Ma 
number of about 7600. Although these experiments utilized a different geometry' Chon was 
used here, and a fluid with a much higher Fr number, they present a useful guide for these 
studies. 

In order to evaluate Ha, che velocity and temperature distributions In the DFZ geometry used 
here were computer modeled with our finite difference techniques, and we found Ma • 6300 for 
this experiment. Previous work with ttoH 03 disks showed this approach is effective /7/> 

It might be noted that one could evaluate Ma more directly from experimental data by moving 
the thermocouple to the center of Che molten disk. This was done on initial experiments, 
and agreement of the Ma Values determined by both techniques was good. In che present 
experiment, che thermocouple located st r ■ 3 mo, because it was felt Chat che amplitude of 
temperature osclllaciona would be higher there chan at che center position. 

CONauSIOHS 

Mathematical models of Marangonl flow in a tin disk heated from below indicaco the flow has 
a two cell structure In which one cell is a chin surface boundary layer which decreases In 
thickness as the parameter (d/t)®Ma/Fr Increases. Finite difference calculations give radial 
distribution of surface velocity which exhibits double peeks as the Marangonl number 
Increases. 

Experiments, In which the surface of the thin nolten tin disk was atomically clean, show no 
evidence of thermal or hydrodynamic oaclllatlona in the melt at Marangonl numbers up to 
6300. The effect of natural convection In these experiments was very small because che 
Bond number, d^gBp/|do/dTj , was kept below approximately 0.05 by maintaining d leas than 
or equal to 0,12 cm. 
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CONVECTIVE DIFFUSION IN ZONE REFINING OF LOW PRANDTL 
NUMBER LIQUID METALS AND SEMICONDUCTORS 
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Several elementary aspects of mass diffusion, heat 
transfer and fluid flow are considered in the context 
of the separation and control of mixtures of liquid 
metals and semiconductors by crystallization and 
float-zone refining. First, the effect of convection 
on mass transfer in several configurations is con- 
sidered from the viewpoint of film theory. Then a 
nonlinear, simplified, model of a low Prandtl number 
floating zone in microgravity is discussed. It is 
shown that the nonlinear Inertia terms of the momentum 
equations play an important role in determining 
surface deflection in thermocapillary flow, and that 
the deflection is small in the case considered, but it 
is Intimately related to the pressure distribution 
which may exist in the zone. However, thermocapillary 
flows may be vigorous and can affect temperature and 
solute distributions profoundly in zone refining, and 
thus they affect the quality of the crystals 
produced. 
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Introduction 


The basic idea of zone refining is that a liquid region or zone is 
created by melting a small fraction of the material in a relatively 
long solid charge, Ingot or feed stock. By moving this liquid zone 
through the charge in one direction the solid phase can be purified 
as the forward surface is melting and the rear one solidifies. 

This is referred to as zone refining. If the liquid zone is passed 
through the charge in both the forward and reverse directions, a 
uniform distribution of Impurities may be obtained. This is zone 
leveling. 

Pfann (1^) first described the essential features of zone 
refining and pointed out its potential as a separation technique. 

In the early 1950* s it was used to provide high purity silicon and 
germanium for semiconductor applications. Since then it has been 
used in a variety of applications. The real power of zone refining 
is that one can pass a molten zone through a solid phase numerous 
times without difficulty. Up to a limit each pass Increases the 
purity of the solid phase by decreasing the concentration of the 
solute. It is by means of multipass operation that great purity can 
be achieved in the solid phase. 

Zone refining is a dynamic nonequilibriura separation process. 
However, to understand its essential features it is important to 
understand the equilibrium concepts on which it is based. Of 
central Importance is the notion of the equilibrium distribution 
coefficient, kg = 03 / 0 ^, where C 3 and are the concentrations of 
solute at equilibrium in the solid and liquid phases respectively. 
Since the value of k^ may be dramatically different from unity, it 
is clear that at equilibrium a solute may distribute Itself between 
solid and liquid phases with a great preference for one or the other 
at a given temperature. On this basis the relative amount of solute 
in each phase can be controlled, and a separation can be carried 
out . 

Phase diagrams, which describe the equilibrium relations that 
exist between phases in mixtures, are often very complex as is the 
case with the nlckel*~alumlnum diagram for example. However, if one 
restricts attention only to a relatively narrow range of concen- 
tration and, in particular, to dilute solutions, then one can 
simplify this description of these relationships markedly by using 
straight line approximations as shown in Figures l(a, b) . 

Figures l(a, b) represent phase relationships for cases in 
which solid solutions exist, a common occurrence in metallic and 
semiconductor systems. In such cases solidification does not cause 
complete separation and the degree of separation depends not only on 
the equilibrium relationships represented by Figures l(a, b) but 
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(a) (b) 


Figure 1* Possible solid-liquid phase diagrams. 

also on the convective-diffusive characteristics of the system. If 
the addition of solute lowers the melting point as in Figure la, 
then Uq < 1; on the other hand if the melting point is raised by 
adding solute, > 1 as in Figure lb. 

A solute distribution exists in the melt because the solidifi- 
cation is carried out at a finite rate. For example, if kQ < 1, 
then solute is rejected and accumulates at the surface which is 
solidifying, and this creates solute gradients in the melt which 
tend to be relaxed by molecular diffusion and any convection which 
may exist. The interfacial distribution coefficient, k, refers to 
the solid to liquid solute concentration ratio at the Interface. It 
is k which is used in transport calculations when one is trying to 
understand the dynamic behavior of zone refining systems. It usu- 
ally is found that equilibrium exists locally at the solid-liquid 
interface, in which case k «= kQ. 

The concentration, C^, of solute in the liquid in the neighbor- 
hood of the solid-liquid interface is strongly influenced by k which 
thereby affects the concentration in the solid phase after solidifi- 
cation in at least two ways. First, the smaller k is the faster 
solute will accumulate on the liquid side of the interface. Second, 
it is the value of Cj_|(0) which determines Cg at the Interface for 
a given value of k. With no convection in the melt, it will be 
shown that the maximum interfacial concentration is approached as 
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the process progresses nnd is, Cj^(O) - Cj_^('»)/kj where C,(«) is the 
concentration of solute In the solution outside the diffusion 
boundary layer. 0^(0) climbs from Cj^(«) to the Cj^(»)/k asymptote. 

The motcrlal to be purified often Is very reactive^ and It Is 
difficult to find container materials which It will not atta>ik. In 
such cases floetlng^zone molting Is attractive and Is used for 
example to grow oxygon free silicon crystals. In this case the 
molten zone Is held in place by Its own surface tension which works 
ag&lnst the action of gravity. A molten zone can be established 
between two rods, one a polycrystalllnc charge or feed rod ond the 
other a crystalline rod of purified material which Is either the 
product of the separation or may serve as the feed for another pass. 
The molten zone may be created In several ways including radio- 
frequency Induction heating as is used for silicon or electron beam 
heating which Is favored for refractory materials. 

The gravitational field limits the length of the molten zone 
that can be sustained by surface tension forces, ond this also 
limits the crystal diameter. The possibility of reducing these 
constraints by decreasing significantly the Intensity of the 
gravitational field has generated considerable interest In 
performing floating zone experiments in an orbit around the Earth. 

Heating from an external source creates rather large tempera- 
ture gradients along the surface of the floating zone which give 
rise to strong thermo-capillary convective flow which mixes the 
melt. Furthermore this flow may couple with natural convection 
flows driven by the gravitational field. These flows give rise to 
time-dependent behavior which is of considerable theoretical and 
practical interest because It creates growth strlatlons In semi- 
conductor crystals which affect their performance, and because its 
origin, and the parameters which characterize it, are not well 
understood. 

A detailed review of zone melting and Its applications has been 
given recently by Shaw (^). In the present paper we shall confine 
our attention primarily To the convective-diffusive characteristics 
of such systems, and we shall strive primarily to obtain a sound 
quali*.ii£ive understanding of their behavior. 

The flow phenomena involved in zone refining will be discussed 
briefly. In particular we shall consider surface tension driven 
flow in a cavity containing a low Prandtl number, Pr, fluid (a low 
Fr number is typical of liquid metals and semiconductors). It will 
be shown that simplified models of such flow, which simulate the 
melt configuration in zone refining, predict multiple steady-state 
solutions to the Navier-Stokes equations exist over a certain range 
of the characteristic parameter. 

The Equations of Motion, Energy and Diffusion for Molten Zones 

In actual zone refining operations one may encounter three- 
dimensional simultaneous heat and mass transfer with moving inter- 
faces, and the system also may be time dependent. The mode of 
heating affects system behavior significantly, as does radiation, 
when materials with high melting points are involved. All of these 
factors complicate the analysis of the diffusional aspects of this 
separation process. Tlierefore we shall treat some rather simplified 
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vcratons of real oyotems In an effort; to gain some Inelght into 
several of the Important phenomena which occur. 

First, VC shall assume the syslem is two-dimensional and that 
the physical properties, except surface tension and density (in the 
body force term of the momentum equations) are constant. With these 
assumptions the continuity equation, momentum equations, energy 
equations and the diffusion equation are given by 
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The boundary conditions for Equations 1-5 depend on the mode of 
heating and the way the separation is conducted. To gain some 
insight into how the separation occurs at the solid-liquid interface 
we first consider Equation 5. We shall assume that the melt is 
quiescent so that the process is governed entirely by diffusion and 
u and V are zero. If there is no convection in the melt and 
directional solidification is occurring only in the x-dlrection, 
then Equation 5 becomes 
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BC 3 C 



( 6 ) 


Let us now assume, as is customary, that the charge is being moved 
through tho heater at a constant velocity U, and that this causes 
the solid-liquid interface to move at the velocity U. Then it is 
convenient to use a coordinate X which moves with Che Interface such 
Chat 


X “ X - Ut 


and Equation 6 becomes 


(7) 
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At time t 0 Che solute concentration in the charge was Cq. 

Since diffusion in the liquid phase is slow, we can assume the dif- 
fusion layer near Che interface is thin compared to the length of 
the melt; and therefore C(C,») " Cq> Now we must consider the 
boundary condition at Che interface* To do this we equate the 
fluxes in the solid and liquid phases at X •* 0 and remember that 
with respect to X there is an apparent convective velocity equal to 
-U* Furthermore, » Ds bo that solid phase diffusion can be 
neglected. Tlierefore 


3C (1-k) U 

(t,0) C(t,0), (9) 

ax 

Cg(t,0> 

where k ■ . If we assume k ■ k„ and use the idealized phase 

CL(t,0) ° 

diagram with straight lines shown In Figure 1, k Is constant and 
this makes it much easier to solve the problem posed in Equation 8 
and 9 together with the other Initial and boundary conditions. 

Order of Magnitude Considerations 

There are various ways of using Equations 8 and 9 to obtain infor- 
mation about the solidification process. The simplest one Is to do 
an order of magnitude analyses, OMA, of these equations. This 
yields immediately that on a relative basis the first, second and 
third terras are of order 1/t, U/fi and Di/6^, where 6 is the approxi- 
mate thickness of the diffusion boundary layer. Equating the first 
and last terms gives 


fi " HoTT, 

L 

and the second and last ^erms give 

S ~ Dj^/U. 

A similar OMA of Equation 9 yields 
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If one combines Equations 10, 11 and 12, one gets 
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Equations 13 and 14 provide a qualitative picture of the relation- 
ship between the concentration at the solid-liquid Interface and 
that In the bulk of the melt. 

We see from Equations 10-14 that two time regimes exist In the 
problem. The first period Is a transient period which exists at 
t « The second period is a steady state one which occurs 

when t » Dl/u 2. Since 1b a physical property which cannot 
be manipulated, we see that the greater the velocity of solid-liquid 
Interface the more qu.^;ckly a steady state Is reached In the system. 
This conclusion is of practical importance because one often wishes 
to distribute the solute throughout the solid with a constant 
concentration. 


Exact Qne-Dlmenslonal Solutions of Diffusion Equation 

In the steady state 9C/ 9t » 0, and Equations 8 and 9 can be solved 
easily to give 


1-k UX 

C /C^ = 1 + exp ( ), 

k D, 


(15) 


which gives the steady state solution for the interfacial concen- 
tration as 


Cj^(O) 
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Equation 16 shows clearly that the solid phase steady concentration 
is the initial melt concentration. It is interesting to note 
that Equations 14 and 16 are identical. 

The solution of Equation 8, if one includes the unsteady state 
effect is more complex, and was apparently first shown by Smith et 
al. (3_) to be, 
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where h ■ (l-'k)U/DL. 

Equation 17, evaluated at X “ 0, together with the definition 
of k, and setting t ■ x/U, enables one to calculate the solid phase 
distribution as 

C_ , / Ux Ux /lix 

- |{l + erf / + (2k-l) exp[-‘k(l-k)— ] erfc[(2k-l) / ]),(18) 
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where x is the distance from the point at which the first solid was 

frozen. 

Effect of Convection on Segregation 

The preceding discussion assumes that no convection exists in the 
melt, and this is rarely, if ever, the case. Next we shall consider 
two approaches which account for convection in the melt, a transport 
mechanism which is especially Important in mass transfer because 
is small and even weak convection markedly alters solute concen- 
tration profiles and may cause macrosegregatlonu First we shall 
discuss film theory which is a very simple approach that gives 
qualitative information and often provides considerable physical 
insight into the >;jachanisms Involved. Second, we shall discuss a 
simplified model of zone refining. 

In film theory one assumes that the concentration changes near 
the interface occur in a very thin region of thickness 6. In this 
region the differential equation which describes the concentration 

is given by the steady state (-g^ 0 ) form of Equation 8. The 

magnitude of 6 is determined by the convection which exists in the 
system. 0n'.‘ estimates 5 by solving simplified convection problems 
such as natural or forced convection along a flat surface, or two or 
three-dimensional stagnation flow, etc. From these simple cases one 
calculates the mass transfer coefHclent, 1%, and 6 is defined as 
fi = essence, one neglects all convection when calculating 

C in the stagnant film, and includes all convective effects in its 
thickness, 6. The stronger the convection the smaller is 6. 

The steady-state solution satisfying Equations 8 and 9 and 
C = Co at X » 6 is 
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Therefore, one obtains the well known Burton-Prim-Slichter (^) 
equation 
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and the problem la reduced to choosing an expression for 6 which 
represents a realistic estimate of the convective patterns which 
exist in Che melt. Let us use a few relevant examples to illustrate 
how film theory may be applied to zone refining. 

By definition, the mass transfer coefficient is given by 


SC(0,y) 
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where y is the distance parallel to the surface and x is iiormal to 
It. Therefore the definition of 5 is, 
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where Sh^ is the Sherwood number defined by ~~ and is 

calculated from a problem similar to the one to which we are 
applying film theory, but one that can be analyzed more easily. For 

example, to calculate one would not Include the moving 


boundary because that is Included in the film model, the steady 
state form of Equation S. One might also use a constant concen- 
tration or constant gradient boundary condition at the Interface 
rather than Equation 9 which also is included in the film model. 

Suppose natural convection is the dominant convection mode in 
the melt. If a natural convection boundary layer created by concen- 
tration gradients exists, then its behavior depends on whether the 
interface Is vertical or horizontal to the earth. If it is verti- 
cal, one can shov' that, 


1/4 

Sh « 0.54 Ra 

y y 


0gACy3 
0.54 ( ) 


( 21 ) 


may be a reasonable approximation, Kays and Crawford (^), and 


5 s 


1 

0.54 


y n y 1/"* 

c— . 

SgAC 


( 22 ) 


If it is horizontal, Stewartson (^) and Gill et al. ( 7 ) showed that 
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1/5 

Shy a 0.75 R3y (23) 

may apply If the interface faces upward and the density of the fluid 
adjacent to it increases with dlatance from the interface; or if it 
faces downward and the density relationship is reversed. In this 
case 

V D. y 

_k_) . (24) 

0gAc 

The important point made by Equations 20, 22 and 24 is that 6 
is a function of y, the distance parallel to the interface, and this 
leads to segregation in the y direction (a nonuniformity of solute 
concentration) in the solid phase as given in Equation 20. 

It is most desirable for 6 to be constant, and there are ways 
to make this happen. If the solid-liquid Interface is circular, 
which is most often the case because the charge is a rod, then 6 can 
be controlled by rotating the rod so that the Interface behaves as a 
rotating disk with angular velocity m. This configuration is widely 
used in analytical and electrochemistry because 6 is essential 
constant if this mechanism controls. Levich (8) has shown that for 
a rotating disk 

D 1/3 

6“ 1.6l( — ) v^v/o). (25) 

V 

Equation 25 shov/s that 6 is constant and its magnitude can be 
controlled by changing w. 

Natural convection to blunt bodies such as cylinders (2- 
dlmensional) and spheres (3-dlmensional) has been studied by Acrivos 
(£) and from his analysis one can show that these configurations are 
characterized by constant boundary-layer thicknesses. For 2- 
dlmensional bodies, 

1/4 

1 R V D 

5 = ( k) , (26) 

0.54 pg4C 

and for 3-dimensional ones, 

1/4 

1 ^ 

« mr ( — ) , (27) 

0.54(2) 8gaC 

where R is the radius of curvature of the surface. Note that 
Equation 26 is identical to 22 except that the length scale is R 
which is a constant, and thus 5 is a constant. This Implies that 
the Interfacial concentration of solute is unifci/m across the 
stagnation region. These results may apply to an interface which 
faces downward (in the direction of the gravity vector) into a fluid 
whose density increases with distance from the Interface or one 
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which faces upward into a fluid of decreasing density with distance 
from the interface. 

When strong temperature gradients exist, natural convection may 
be primarily induced thermally or both heat and mass transfer may 
play comparable roles. In these cases the situation is more 
complex, because the n'jinber of parameters increases. In liquid 
metals and semiconductors the Schmidt number, v/D^, is several 
orders of magnitude greater than the Prandtl number, v/a, and this 
enables one to solve for the concentration profile in a rather 
general way without great difficulty as will be discussed next. 

The number of alternative configurations of the melt and modes 
of heating it that may exist in zone refining is very large. There- 
fore it seems desirable to have a method which enables one to esti- 
mate macrosegregation under a wide variety of flow conditions. Such 
an approach is not available now, but some progress toward it can be 
made by noting that liquid phase diffusion is characterized by a 
large Schmidt number, which implies that diffusion boundary layers 
are thin compared to momentum boundary layers. It seems that 
Llghthlll (1^) was the first to suggest that one can restrict 
attention to the velocity field near the interface under these con- 
ditions, and by doing this one can derive the expression 


0.894 (9D 
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1/2 , 1/3 

dy] . 


(28) 


Equation 28 is a general result for high Sc, two-dimensional flows 
in which the diffusion boundary layer thickness is zero at y *» 0. 

It includes Equations 22 and 24 as special cases, but it does not 
apply to systems in which 6 0 at y = 0, such as stagnation 

regions, and also it does not Include Equation 25. Equation 28 can 
be applied to flows which are driven by temperature differences 
regardless of the magnitude of the Prandtl number. 

To use Equation 28 for thermally driven free convection 

boundary layers one simply calculates « -p from known 

solutions. Then this result is inserted in Equation 28 and 6 in 
Equation 19 or 20. 

Unfortunately, because of the variety of factors, such as 
shape, mode of heating and orientation, that is possible in melts, 
and the complexity of the flow patterns which may exist, it is 
extremely difficult to offer general rules, a priori, on how to 
estimate 6. One needs to examine each particular case carefully to 
obtain even a qualitative understanding of the macrosegregation that 
may occur in the crystals being produced. However, one procedure 
which seems to yield generally beneficial results is crystal rota- 
tion as predicted by Equation 25. 


Simplified Model of Surface Tension Driven Flow in a Two-dimensional 
Molten Zone Supported on the Bottom 


The preceding discussion showed Chat steady-state natural convection 
often leads Co undesirable macrosegregation. It also has been shown 
by Gill (11) that natural convection flows may become unstable, and 
Carruthers (12 ) and Milson and Pamplln (13) have discussed the 
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Implications oil the resulting oscillations on crystal growcht In 
this section we shall examine exact solutions of the Navler->Stokes 
equations for a two-dimensional simplified model of a molten zone 
which is in the form of a cavity or slot of llquld^of depth d, 
supported on the bottom, but with a free surface on top> The zone 
is heated over the length, -il < x < i, by a flux, q, and cooled on 
its ends at x » IL, where L > We shall study the core region 
inside -i < x < i for which a similarity solution exists. Thus we 
arc neglecting end-effects. 

One can show that Equation k is satisfied by a temperature 
distribution in the form 

0 „ 1 — ^ • - gi(n) + 82 (n)(-) , (29) 

’^hot “ ^cold ^ 

where and g 2 are functions which are determined from 

g' ' + 2A^Ma f g' - -2A g (29a) 

I 1 2 

and 

g" + 2A3Ma[f g' - f g ] - 0 (29b) 

2 2 2 

with the initial conditions 

g'(0) - g’(0) - 0 (29c) 

1 2 

8 (0) » -8 (0) - 1, (29d) 

< 1 2 

and T) = y/d. Here f( n) and f'(n) are functions related to u and v 
which will be determined later; Ma is the Marangonl number 
arilda/dTl d 

; and A is the aspect ratio — . In the limit Pr + 0, the 

po il 

solution given by solving Equations 29 implies a constant heat flux 
along the bottom of the cavity, n " 1» and a zero heat flux from the 
free surface, n ” 0, into the vacuum surrounding the liquid zone. 

For non-zero Pr the flux at n = 1 varies with x. Obviously Equation 
29 implies that the liquid surface temperature varies as 

X 2 2 

t(x, 0) = Tj^- (Tj^-tJ H = T^- ATx , (29e) 

where T^ is the temperature at x = 0 and T^, is that at x “ ±2.. 

Gill et al. (1^) have shown by numerical computation that Equation 
29e is a good approximation to a constant heat flux for fluids with 
finite values of Pr which are typical of liquid metals. Therefore 
the following discussion applies best to low Pr fluids, a category 
that Includes essentially all fluids of Interest in semiconductor 
technology as well as all liquid metals. On the other hand, the 
most complete data on thermocapillary flows in molten zones has been 
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reported by Frelsser^ Scharmann and Schwabe (15) and Schwabe and 
Scharmann (16) ‘ for NaNO 3 which has a Prandtl number o£ 9 and for 
heating from the side which Is not consistent with the present 
theory . 

Levlch ( 8 ) has discussed capillary motion In two-dimensional 
creeping flows In which the surface was flat> Ylh ( 17 ) pointed out 
Inconsistencies In Levlch* s analysis which were associated with the 
assumptions of a linear distribution of surface tension with 
distance along the Interface, and with the deflection of the surface 
which Inevitably occurs when capillary flow exists. He noted that 
under certain circumstances steady flows may not exist. Ostrach 
(18, 19) has discussed scaling problems in capillary flows. 

Recently, Sen and ^avls (20) studied capillary flow in bounded cavi- 
ties in which I Is small, end effects are present, and the flow Is 
very slow and the cavity Is heated from the side. Cowley and Davis 
(21) studied the high Harangonl number "Thermocapillary analogue" of 
a buoyancy driven convection problem solved by Roberts (22). Later, 
we shall make some comparisons between our results for the 
deflection of the surface and those of Sen and Davis (20). 

The boundary conditions for Equations 1-3 which will be satis- 
fied at X ■ 0 and the solid bottom are straightforward, but those 
for the free surface are rather complex and call for some dis- 
cussion. At X » 0 we have a stagnation condition and at y » d there 
is no slip and no penetration of fluid through the solid bottom. 
Therefore 

u(x, d) « v(x, d) = 0, (0<x<J,) and u(0, y) = 0, (0<y<d) (30) 

The free surface Is not flat in general and the boundary con- 
ditions on this surface require careful consideration. The kine- 
matic condition at y -h(x), where h(x) is the deflection from the 
flat surface at y = 0, is 


dh 

v(x, -h) u(x, -h). (31a) 

dx 

Equation 31a is the steady state form of the kinematic condition 
which also Is used to describe wave motion as discussed on page 595 
of Levlch* s book, Physicochemical Hydrodynamics . One also must 
equate the normal and tangential components of the forces in each 
phase at the free surface. Since we consider a gas-liquid inter- 
face, we neglect gas phase resistance due to Its viscosity and 
Include only the pressure it imposes on the Interface on the gas 
side. Therefore, at y = -h(x) the tangential component of the 
stress tensor for the liquid phase is equal to the tangential force 
created by the change In surface tension with temperature In the x 
direction. Thus the tangential force balance at the interface 
becomes 


-y 

1 + h2 
X 


2 3u 3v 3v 3u 

[(l_h )(_ + _) -2h^( )] = 

^ ^ ax 9y ax (1 + h2) 

X. 


(31b) 


and the normal force balance becomes 


XX 


2 y 3v 3u 3u 3v oh. 

- (P-P„ ) + [— + h — + "I” 

® (i + h2) ay ’‘ay * ax ax (i + h^) 


T77* (31c) 


where a is the surface tension and the x-subscripts on h and c 
denote differentiation with respect to x, and pg Is the pressure 
In the gas phase over the liquid which at zero g is taken to be 
zero. 

To solve Equations 1-3 together with 30 and 31a-c we take note 
of the important role played by the surface tension number, e ■ 

where Oq Is the surface tension at temperature Tq. In the 
liquid metal and semiconductor systems we are considering e Is very 
small (i.e. e « 10""^/d for liquid tin); therefore a first order 
perturbation solution In e should be adequate. Consequently we 
introduce a stream function, i|< and a similarity coordinate, q, so 
that Equations 1-3, 30 and 31a-c are satisfied to 0(e), and the 
computations are reduced to solving a two-polnt-boundary-value 
problem for a single nonlinear ordinary differential equation in the 
dimensionless stream function f. The procedure for finding the 
base flow solution, which is equivalent to finding f and q for e = 

0, is similar to that employed by Carter and Gill (23 ) to treat the 
entirely different physical problem of flow with heat transfer, 
blowing or suction, and natural convection in flows between 
parallel-plates or in tubes. As usual we define 


a4< 34' 

u = , “V “ • (32) 

3y 3x 

We choose q « y/d which suggests 

q u(x,y) 

4' = Up(x) d / dq. 

^ o U^(x) 

If u(x,y), scaled by Ug(x), is a function of q only, then 


4- = Uj^(x) d f(q). 


(33) 


Equations 2, 32 and 33 lead to 
d2 dU 

f. M + « (ffll _ ft 2) 

V dx 


which satisfies Equations 1-3 exactly. Clearly the coefficient — 
dUn 

— ^ is constant, and Equation 31b is satisfied, if 
dx 


d2 

3p 



— “ 6 

(34) 


dx 

d2 
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(34a) 


so that 


u - 2 (v/d) A^Re Xj f (n), V - "•2(v/d) f (n)* 35(a,b) 


where x^ >■ x/ £ a^d A > d/ Z> To obtain Equation 34a we use Equation 
29 and o » + ^(t-Tq ) which yields a/ o^ ■ 1 + ARec x^. Equation 


34 then becomes 


d 3 Ma 

f" ' + 2 (-) — (ff ' - f 2) 
£ Pr 


0 , 


(36) 


and the boundary conditions become 


f [O ) " f (l ) " f ’ (l ) ” f ' ' (o ) + 1 “ 0, 


(37) 


where 


Ma 


Pr 


Re 


AT£ |da/dTl . a , 

( j and the last boundary condition 

ya V 


given by Equation 37 is obtained by substituting Equations 32a, b 
into Equation 31b which is satisfied to 0(e). Since Equation 36 is 
of third order and 37 includes four boundary conditions, the 

constant 6 can be determined as a function of Q » 2(-^)3 and the 

solution for f can be related to u and v via Equations 32 and 33. 

An important feature of Equation 36 is that it retains the essential 
nonlinear features of the Navler-Stokes equations. It will be seen 
later that the nonlinear terms are of great importance even at rela- 
tively small values of Q, and it is they which determine if multiple 
solutions to Equation 36 do exist, 

Brady and Acrivos (24) studied an entirely different physical 
problem with different boundary conditions from those given in 
Equation 37 and cautioned that the numerical solutions they obtained 
for a system with a capped end differed from the similarity 
solutions. Their problem differs fundamentally from the present one 
in several ways. First, the flow in our problem is symmetric about 
X 0, because surface heating is symmetric about x » 0 over the 
region -£ < x < £, and this causes u = 0 at x * 0. Hence, our 
solution is exact near x 0. Their problem is antisymmetric about 
X e= 0, and they apparently have a "collision region" near x = 0 
which affects the structure of the entire flow. Second, Brady and 
Acrivos specify the shape of their cavity, h(x), as part of the 
problem statement, and their choice of shapes may be inconsistent 
with the pressure distribution in their similarity solutions which 




they use as the criterion for agreement with their numerical calcu** 
latlons. We solve for the shape of the zone and require It to be 
consistent with the pressure distribution to 0(c) In Equation 31c> 
Third, the shape of the zone affects the boundary conditions on v 
through the kinematic condition, and therefore our boundary con- 
ditions are different from theirs* Furthermore, Gill et al> (1^ ) 
found good agreement between similar solutions for thermocapillary 
flow In an Infinite disk and numerical solutions to a finite dl')k 
with end effects. In the thcrmocaplllary problem (1^) the boundary 
conditions were Identical to those In Equation 37. 

We studied the two- point-boundary-value problem (TPBVP) 
Equations 36 and 37 (corresponding to a slot with one free surface) 
by suspending the trivial differential equation 

B' “ 0 (38) 

to 36 and using the TPBVP solver BOUNDS, P. Deuflhard (^), to solve 
the resultant well-posed problem, Equations 36, 37 and 38. BOUNDS 
Is an accurate, multiple-shooting code. We began with the solution 

f = n(l-n)2/4, B - 3/2, (39) 

which corresponds to Q » 0. We found a family of solutions for 
Q > 0 evolving from Equation 39; we first slowly Increased Q and 
then increased Q In larger increments up to approximately Q - 8000. 
For continuation we used linear Interpolation of the solutions 
obtained for Qj and Qg (Q^<Q 2 ) provide the guess for the solution 
at Q « Q3>Q2' Above Q » 8000, B began to change more rapidly, and 
we replaced (38) by 


Q' “ 0, (40) 

In order to prescribe B Q^d compute Q and f. We continued by linear 
interpolation (as with Q) until B “ -0.75 where Q " 20,000. A 
summary of the numerical results Is provided In Figure 2. 

For Q 1 6900 the solutions had 2 cells; see Figure 2. Above Q 
" 6900, a third cell appeared and Increased In strength as B was 
decreased; see Figure 2. A local maximum In Q was reached at Q » 
8614.48, B = -0.0675; a local minimum in Q was reached at Q “ 
8128.32, B = -0.21. Thus for 8128.32<Q<8614.48 we found three dis- 
tinct solutions of the TPBVP, Equations 36 and 37. Figure 3 com- 
pares the three solutions of Q « 8600 with that for Q “ 500. It is 
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Figure 2. Results of numerical evaluation of TPBVP, Eqna. 36 and 
37. 



Figure 3. Numerical predictions of axial velocity versus height 
above the support surface. 
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seen that the dlfferencco between the velocity distributions on the 
middle and lower parts of the branch arc so small as to be indls** 
tlngulshnblc on the figure. 

The rapid change in 3 with Q for 0<Q<10^ shows clearly that 
nonlinear effects are Important in Equation 36 even when the 
Marangoni number is relatively small. 

Now let us calculate the deflection of the surface and 
determine how it depends on the parameters of the system. To do 
this, substitute Equations 35a, b into 31'c, together with 
h •• h^ + hjC, and by equating coefficients of e one obtains hQ ■ 

0 and for - h j/d the equation 

d2n, p d2 3 2 

— ■ [4ARe f’ (0) + + — Re x j] (41) 

dx ^ UV A 

with the conditions that the deflection at the ends is zero and that 
the liquid Just fills the cavity. To avoid the complications 
associated with the end-effects we take i •• L for the purpose of 
calculating H Thus 

i 

Mi(-a) - Hi(i) - / d Xj - 0. (42) 

-I 

Then the final result for h is 

3 Re 

h/d e(1 - 5x,2)(i - X 2), (43) 

60 A 

Equation 43 shows that the deflection is small because e is 
small, that it increases with Re/A and that its sign depends on the 
sign of 3. However, it must exist in order to make the pressure 

d 2 8p 

distribution given by =3 consistent with the boundary 

MOo 8x 

condition given by Equation 31c. 

Sen and Davis (20) have performed an interesting analysis of 
capillary flow in a slot heated from the side (a half zone) by using 
the matching procedure suggested by Cormack, Leal and Imburger (26) 
for natural convection in a slot. In this way they take end-effects 
into account. However, they assume that the nonlinear terms in the 
momentum equation can be neglected and perform an analysis for small 
values of d/i, which they use as their perturbation parameter. 

Their result, to the first order, is 

h Re 

— “ e[x, - 3 x.2 + 2x,2], (44) 

d 8A 

where x^ Is Che dimensionless distance from the hot wall in their 
case. 

One cannot compare Equations 43 and 44 quantitatively because 
they apply to different configurations and therefore Che distri- 
bution of h with X is different. However, it is interesting to 
compare these results qualitatively. Let us look at Che main point 
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of similarity first* In tho limit Ro 0| 0 3/2, and for this 

value of 6, Equation 43 predicts a maximum deflection of 0*023 Rec/A 
while Equation 44 predicts 0*012 Roc/A* Obviously, tho parameter, 

I’afl ^ 

Ree/A ■ — r '7 - on which the deflection depends la the some In both 

Oo(d/i)Z 

cases. We note also that this parameter Is Inversely proportional 
to (d/i)^' Thus the smaller (d/Jl) is, the larger the fractional 
deflection, h/d* On the other hand In the linear flow case Equation 
43 predicts a larger deflection than Equation 44; this may be due to 
the neglect of end^cf facts In our case or the fact that heating from 
the side Is studied by Sen and Davis and we have considered surface 
heating. The most Important difference between Equations 43 and 44 
is the appearance of 0, which depends on the parameter 2(d/P.)^ Re 
and therefore reflects the nonlinear terms In the Navler-Stokes 
equations. For a fixed (d/l) one secs from Figure 2 that h/d first 
decreases with Re because 0 decreases rapidly, and then h/d 

Increases with Re because 0 remains relatively constant up to Re 

-6000. 

Figure 2 shows that valuet of 0 change rapidly from positive to 
negative as 2 (^) Re increases from zero to a little over 1100. 

Beyond a value of 2('^)\e 6900 a third cell appears, and at about 

8100 multiple solutions for 0 appear after :it has changed very 
slowly from ab(^ut 1100 to 8000. The rapid changes in 0 at small 

values of 2('^) Re suggest that nonlinear effects play a decisive 

role In the behavior of the deflection, h/d, in Maratigonl flow. 

We speculate that the existence of multiple solutions to the 
Navier-Stokes equiitions for surface tension driven flow in a cavity 
is evidence of Ijistabllity in this system, and this suggests the 
need for a stability analysis of the velocity profiles. It should 
be noted that the simplified model of the molten zone which we are 
using does not Include the blocking effect of the solid-liquid 
Interface at x - ±1. 

Calculations by Gill et al* (14) , in which a similar simplified 
model was used for a disk shaped molten zone, agreed well with more 
detailed finite difference calculations, except in the immediate 
neighborhood of the solid liquid Interface* Both types of calcu- 
lations for the molten disk suggested that only one steady-state 

A 3 

similarity solution exists and that 0 changes rapidly with 2(^) Re 

so that nonlinear effects also are important in chat configuration. 

The disk problem was reduced to equations similar to Equations 
36 and 37 except, for the disk, the term ff*' is replaced by 2ff " . 

Therefore, the term v^ is weighted more heavily compared in the 

disk than in the cavity and this eliminates multiple solutions* 

This seems to mean that a relative increase in the convection of x- 
direction momentum normal to the free surface compared to its con- 
vection parallel to the free surface favors stability in the zone. 
However, a stability analysis is required to confirm this hypothe- 
sis. 
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Smith and Davis (^) have carried out an inccreacing stability 
analysia on liquid loyers in a slot heated from the side. They 
ncBicctcd end effects and the nonlinear Inertia terms in their basc~ 
flow velocity distribution and identified three types of instabili- 
ties: stationary roll cells, hydrothermal waves and surface waves. 

In systems with return flows, as is the case in floating or 
supported zones, stationary roll cells were shown not to exist in 
slots heated from the side. Surface waves were found to be the 
prevailing cause of instability at very low Fr. However, it seems 
tenuous to try to extrapolate their results to the return flows with 
significant inertial effects and surface heating which we studied 
because the nonlinear inertia terms have a profound influence on the 
velocity profile, the surface deflection and the pressure distri- 
bution; therefore they would influence strongly all modes of insta- 
bility. 

It is cleat that vigorous thermocaplllary flows with surface 
velocities on the order of 10 cm/scc nay exist in low Frandtl number 
liquid metal or semiconductor zones. These flows no doubt have a 
profound effect on the temperature and solution distribution in 
these zones, and thereby they affect crystal quality. This con- 
clusion is consistent with that drawn by Wilcox and coworkers, (28) 
and (29), who first called attention to Che potential importance of 
surface tension driven flows in float zone solidification. Schwabe 
and Scharmann (16) have recently summarized experimental investi- 
gations on a wide variety of high Pr number fluids. Surface veloci- 
ties have been directly measured in transparent fluids with Pr 
numbers varying from 0.3 for molten KCl to over 1000 for silicon 
oils and molten glass. The study on molten glass, for example, 
recently reported by McNeil et al. (30) with both top and bottom 
heating of a liquid bridge reported flows of roughly equal velocity 
in the range of 0.1 to 1 cm/s. This study and the several others 
listed by Schwabe and Scharmann (16) confirm the predominance of 
surface-tension flow over natural convection even at 1 g, and it 
seems reasonable that this same conclusion will apply to liquid 
metals and semiconductors at Fr numbers which fall into the .01 to 
.03 range. 
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III. Summaries of experimental work carried out during the final 
reportiiig period, October 1, 1983 to June 30, 1984. 


Experiment 1, Detection of Oscillatory Flow 
Experimental Design 

The aim Af initial experiments was to detect oscillatory Marangoni 
flow In the tin disk floating zone geometry (DFZ) utilizing a thermo- 
couple to detect oscillations produced by oscillatory fluid motion. 
Theoretical models predict values of either a surface Reynolds number 
or a Marangoni number at which one expects to see the onset of oscilla- 
tory flow. Therefore, an Initial phase of our experiments was to util- 
ize the computer program developed by C. C. Hsu (which utilizes finite 
difference techniques to model our DFZ geometry) in order to design an 
experiment which would allow us to achieve high values of (Ma)g, the 
Marangoni number utilized by Schwabe, and (Re)^^^, a surface Reynolds 
number utilized in our theoretical analyses. 

The geometric variables of the DFZ geometery are illustrated 
schematically in Figure 1. Four parameters enter into the design, the 
total tin disk radius, R, the molten tin disk radius, H, the heater 
radius, a, and the disk height, d. The two dimensionless parameters 
are related to controllable experimental parameters as follows: 

|do/dT| 'AT-f 

(Ma)g 

|ia 

where da/ dT « temperature derivative of surface tension, AT “ 
temperature difference between hot and cold sinks, I = radius of the 
molten zone, p •= viscosity and a » thermal diffusivity. 
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Figure 1. The disk float zone (DFZ) geometry. 
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where p •= density, and is a reference velocity defined as, 
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where ATn is a reference temperature difference defined as 
R 
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.kjd 


( 4 ) 


where kjj^ = liquid thermal conductivity and P is the power into the 
disk, which may be evaluated from the heat flow in the solid portion of 
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P 


27f k d (t - T ) 
8 mp c ^ 

Ln(R/ll) 


(5) 


where kg ■ solid thermal conductivity, Tjnp “ melting point temperature, 
T(. * cold sink temperature. 

There are two major limitations upon our experimental system, the 
temperature of the cold sink cannot go below 40 to 50°C and still main- 
tain good control and the high temperature cannot go above 500 to 550°C 
because of limitations Involving a solder Joint in the heater design, 
plus a solder used to attach the heating element. Our design objective 
Is to achieve conditions which correspond to a value of (Ma)g greater 
than 7400 within these limitations on temperature. The value of 7400 
is that value found by Schwabe^ to cause oscillatory flow In NaN 03 
under cylindrical geometric conditions. Based on Initial studies 
(Annual Report 1983) we chose for the optimum conditions within the 
power limitations of our present heater design the following con- 
ditions; a = 3 mm, d = 1 mm, s » 6 mm, and R « 15 mm. The computer 
output for these conditions gave a dimensionless temperature at the 
heater center, r = 0 and the heater interface, Z =■ 0, [0(0,0)] as a 
function of (Re)j^ from which we obtained the temperature at the heater 
Interface, 

= 232 + 6(0,0) . ATj^ (6) 


where ATj^ is given by Eqn. 4, The value of (Ma)g is related to 


(Re)j^ as, 


(Ma)^ - (Re)^ Pr (|) 8(0,l) 


(7) 


where 6(0,1) Is the dimensionless temperature at r ■ 0 and the upper 
disk surface. By combining equations 2 to 5 one obtains an expression 
for as a function of the controlled parameters of the experi- 

ment . 

2pd Ida/dTl k (t - T 
g ^ ' 8 ^ mp c 

^ W^k^ln (r/£) 

From Eqns. 4 through 8 and the computer output for 0(0,0) and 
6(0,1), the significant parameters are given In Table I for the 
conditions, R « 15 mm, £ » 6 mm, d « I mm, p ® 7 gras/cm^, |da/dl| « 
0.17 dynes/cm-^C, p = 0.018 dyne-s/cm^, « 0.305 j/cm-C, " 0.06 
j/cm-C, a = 3.0 mm and Pr « 0.014. 


Table I. Computer Model Predictions 


CnOi, 

GO 

Hot Temp., Tj^ 

Cold Temp, T^, 

Power, P 

10“* 

482 

249 "C 

225“C 

2.65 watts 

5 X lO** 

2200 

314“C 

200"C 

13.2 watts 

10^ 

4370 

388“C 

lei^c 

26.5 watts 

5 X 10^ 

17800 

882‘’C 

- 90 “C 

133 watts 


These values along with some others have been plotted as Figure 2. 
From these plots it is seen that a (Ma)g value of 7400 is achieved at 
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Figure 2. Plots of the data shown in Table I, for Sn with R « 15 mm, 
I 6 mm, d = 1 mm and a « 3 mm. 


a Th value of 480°C and a T(. value of 115“C. There are two experi- 
mental parameters which we can easily control and measure, T^. the cold 
sink temperature using thermocouple number 2 (t.c. //2) on Fig. 1 and 
A, the radius of the molten Sn zone. T^ is controlled by the tempera- 
ture of the hot oil bath. For a given Tj. the value of I is adjusted by 
power input to the heater and it is monitored by means of a telescope. 
Therefore, from the plot of Fig. 2 we may evaluate both and 

(Ma)g for any given value of Tg and the power input conditions adjusted 
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to give H ■ 6 mm. To achieve high values of (R®)j^ (Ma)g one needs 
to lower Tg. Selecting a minimum achievable !'(. of 70®C gives a maximum 
obtainable (Re)j^^ value for this design of 2.5 x 10® which corresponds 
to (Ma)g ■ 10,000 which is just slightly above the most recently 
reported values of Schwabe^ for his NaNO^ fluid. 

Experimental Result s 

Experiments were carried out in our UHV system. Utilizing a cold 
sink temperature of we were able to produce an 11.7 mm diameter 

zone with a clean upper surface at a value of Ma » 4300 (see section II 
for details). This was the maximum Ma possible with our heater design 
(see 1983 annual report for heater design details). There are two 
points in the design, a solder joint and a current limitation on a 
vacuum feedthrough, which set this limit. A new design has been 
developed which should allow Ma value of 10 ** to be achieved. As may be 
seen from Fig. 2 the heater must be capable of raising the molten tin 
temperature to 570“C at the heater Interface, which means the heater 
Itself must operate somewhere above 600®C. 

Our experiments have not detected any temperature oscillations due 
to Marangoni convection. This could mean either (1) temperature oscil- 
lations are in fact not present, or (2) temperature oscillations are 
present, but they occur at frequencies too high for our measuring 
system to detect. Clearly, it was very important for us to evaluate 
the frequency response of our thermal measuring system, and experiments 
were carried out to this effect. 

A frequency generator was set up, attached to our system and 
experiments were carried out under several configurations. It was 


RECIPROCAL OF FRACTIONAL 
ATTENUATION OF AMPLITUDE 
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quickly determined that the best frequency response was obtained by 
feeding the thermocouple output directly into the chart recorder. 
(Initially we had utilized a Kelthly 181 nanovoltmetcr and although the 
output ot this device allowed improved sensitivity it greatly reduced 
frequency response) . Measurements were made on the chart recorder 
using the same sensitivity scale as used for recording temperature in 
the experiments. Sine wave voltages of amplitudes varying from 10 to 
40 \N were measured as a function of frequency from 0.02 to 20 Hz. The 
attenuation of the amplitude was determined and the results are pre- 
sented as the reciprocal of the fractional attenuation of amplitude 
versus frequency on curve A of Fig. 3. 



Figure 3. Response of thermocouple output as a function of signal 


frequency. 


MINIMUM DETECTABLE 
TEMPERATURE AMPLITUDE 
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In gonoral our experimental work sliowa that we can detect a mlnl~ 
mum amplitude of oscillations on the order of 1 pV with our type K 
thermocouples. The minimum detectable temperature amplitude given as 
the right hand scale of Fig. 3 is for this 1 pV detection limit on the 
amplitude of the oscillation. Referring to curve A, Fig. 3» this means 
that at 10 Hz in order to detect a signal oscillation its amplitude 
would have to be > 13.5 pV> which corresponds to an oscillating temper- 
ature amplitude of 0.34'’C with our type K thermocouples, as is shown on 
the right hand ordinate of the plot of Fig. 3. 

This approach assumes that the thermocouple itself has an infi- 
nitely fast time response and we are limited only by the response of 
our recorder system. Experiments on the time response of small thermo- 
couples in molten metal environments are available in the literature^ 
and they show that a 10 rail (250 pm) thermocouple has a time constant, 
T, of 35 ms in molten Sn-Pb. 

This finite time response will also reduce the amplitude of an 
oscillating signal and it may be shown that the reciprocal of the 
fractional attenuation of the amplitude of the signal is given as 
[1 + ( (jJt) 2] ^ where mis frequency in radians per second. This 
function is shown as curve B on Fig. 3. The total attenuation of the 
signal is then given by curve C which is the product of A and B. 

These results predict the following picture. With a 10 mil 
thermocouple having a 35 mS time constant attached to our chart 
recorder system we should be able to detect temperature oscillations 
greater than or equal to the values given by the right hand ordinate 
using curve C. The frequency dependence of these results are summa- 
rized in Table II. From these results it is clear that temperature 
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Table II. Tabulation of results of Figure I, curve G. 


Frequency 

Minimum Detectable Temperature 
Oscillation Amplitude 

<0,1 Hz 

.025*G 

1 

Hz 

.031“C 

2 

Hz 

.04S”C 

5 

Hz 

.155“C 

10 

Hz 

.950“C 


OBcllIatlons above 5 to 10 Hz are unlikely to be detected by our 
present system unless their temperature amplitude is quite large. The 
results also show that to improve the frequency of our system we can 
(1) reduce the diameter of the thermocouple in order to reduce its time 
constant, (2) increase the Seebeck coefficient of the thermocouple, and 
(3) Improve the frequency response of out recorder system. This latter 
improvement could be achieved by going to a digital data acquisition 
system. 

These results have been presented briefly in our COSPAR paper, see 
pages 15-16. We conclude that at Ma =» 4300 in the disk system no 
temperature oscillations occur of frequencies less than 5 to 10 Hz. 


Experiment 2, Interface Shape Evaluation 

The purpose of this experiment was to detect Marangoni flow in 
this opaque metallic system by observing its effect upon the shape of 
the solid-liquid Interface. The procedure would involve a set of com- 
parison experiments. In the reference experiment the top molten Sn 
surface would be cleaned of oxide film with argon ion sputtering, and 
the shape of the solid-liquid interface formed at the point of furthest 
melt-out would be determined. The predicted Ma surface flow velocity 
toward the solid is around 25 cm/s for the Ma * 4300 experiment. Since 
the bottom surface of the disk has a thick oxide layer upon it we 
expect no Ma flow there. Hence, the rapid Ma flow on the top surface 
should cause the S/ £ interface to be melted outward toward the cold 
sink successively further as one goes from bottom to top of the disk. 
Therefore in the subsequent comparison experiments one selectively 
oxidizes the upper surface to different oxide layer thicknesses and 
evaluates If some minimum oxide layer thickness is adequate to maintain 
the S/ £ interface a constant distant from the cold sink from top to 
I bottom of the disk. 

The technique for marking the S/£ interface was simply to use a 
eutectic alloy, Sn-Au. This alloy was chosen because (I) Au is a noble 
metal and might not change the oxidation characteristics of molten Sn 
significantly, and (2) a eutectic occurs at only 10 wt.% Au and it has 
a large volume fraction of second phase, ~ 30%, while only lowering the 
melting point of Sn by IS^C. The microstructure of a eutectic la very 
sensitive to solidification conditions. Therefore, by starting with a 
disk of a given microstructure and melting out to radius R, one expects 
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to be able to easily detect the mclt**out Interface shape by doing 
mctallogrophlc analysis on sectioned disks* 

Techniques were developed to prepare disks of Sn *■ 10 v/o Au 
alloys and two experiments were run with essentially Identical results. 
It was found that after melting out the zone to a given diameter upon 
holding for 10 minutes or so a ridge of solid material began to be 
deposited all along the edge of the solid-liquid Interface. In a few 
hours time the ridge became quite large as Is shown In the sequence of 
photos presented in Fig. 4. Sectioning of the sample revealed that the 
S/1 interface was clearly delineated, but the presence of the ridge of 
solid material was too disruptive to allow meaningful Intcrpt-^tatlon of 
these experiments. If we could have sollii't.*;flcd the samples before 
enough time passed to build up the ridge material Chen the experiments 
would have been pursued. Unfortunately, we could not clean the molten 
surface by sputter cleaning in less than an hour or so. This suggests 
cleaning the surface prior to melting. We found that it was necessary 
to clean the surface after additional solid was melted because oxygen 
and sulfur were detected at the newly melted surface, due apparently to 
liberated atoms from the melting. Hence, Che surface could not be 
cleaned prior to melting and thereby be assured of no oxidation. 
Consequently these experiments were abandoned. The cause of the ridge 
buildup will not be addressed In this report. It does not occur with 
pure tin. 
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Experiment 3, Macroaegregation Study 

The purpose of thia experiment was to detect Marangoni flow in the 
DfZ apparatus by evaluating its effect upon macrosegregatlont As with 
expQriraent 2 this technique involves a comparison study between a disk 
solidified with no oxide and one solidified with oxide deliberately 
controlled to a desired thickness level. The control experiment Is the 
oxide free experiment where one expects very large convection veloci- 
ties along the top surface, V s 25 cm/s with Ma 4300. In this case 
one expects an effective distribution coefficient k,, nearly equal to 
the equilibrium coefficient ko» Py repeating the experiment with an 
r'<ide film of sufficient strength to suppress Marangoni convection one 
expects the value of kg to rise to 1 after an initial transient whose 
length will depend upon the rate of solidification. Hence, by 
measuring kg one can indirectly evaluate the effect of convection and 
hence the effectiveness of oxide films upon suppression of this con- 
vection. In these experiments it was planned to measure kg by util- 
izing the radioisotope Bi^^^ at levels of around 200 ppm to determine 
composition variations as a function of fraction solidified. 

Wa had difficulty procuring some Bi^O? ^ut a source was located 
and an experimental technique was developed to plate out the Bi^O? and 
to dope it into Sn disks at around the 200 ppm level. This work relied 
in part on pur previous studies with this Sn-Bi^O'^ system. In 
addition, it was necessary to be able to control the power input to the 
heater in a programed way in order to solidify the disk at a constant 
rate. We were able to borrow a prograinable power supply for this work 
and a relatively inexpensive Apple lie computer was purchased for the 
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programming. This work, coupled with the fact that our Auger CMA had 
to be returned to the manufacturer for maintenance <2 month turnaround) 
was a sufficient delay that the first Bi^O^ experiment was not run 
prior to our June 30 termination date. 


55 


References 

1. F. Preisser, D. Schwabe, and A. Scharraann, J. Fluid Mech. 126 , 545 
(1983). 

2. D. Schwabe and A. Scharmann, Paper G.1.2.4 COSPAR Conference, 

Graz, Austria, July 1984. 

3. J. Nanlgan, Namao Temperature Handbook p. L27, Namac Corp., 
Framingharo, MA (1981/82). 

4. R. P. Benedict, Fundamentals of Temperature, Pressure and Flow 
Measurement , p. 270, John Wiley, NY (1977). 

5. J. D. Verhoeven, E. D. Gibson, and R. I. Griffith, Met. Trans. 

475 (1975). 

6. J. D. Verhoeven, E. D. Gibson and B. Beardsley, Met. Trans. 

349 (1975). 


6U.S. GOVERNMENT PRINTING OFFICE 1985-644-048/10083 


